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Abstract 

We take an Effective Field Theory (EFT) approach to unifying existing proposals for the 
origin of cosmic acceleration and its connection to cosmological observations. Building on earlier 
work where EFT methods were used with observations to constrain the background evolution, 
we extend this program to the level of the EFT of the cosmological perturbations - following 
the example from the EFT of Inflation. Within this framework, we construct the general theory 
around an assumed background which will typically be chosen to mimic ACDM, and identify 
the parameters of interest for constraining dark energy and modified gravity models with ob- 
servations. We discuss the similarities to the EFT of Inflation, but we also identify a number of 
subtleties including the relationship between the scalar perturbations and the Goldstone boson 
of the spontaneously broken time translations. We present formulae that relate the parameters 
of the fundamental Lagrangian to the speed of sound, anisotropic shear stress, effective New- 
tonian constant, and Caldwell's w parameter emphasizing the connection to observations. It is 
anticipated that this framework will be of use in constraining individual models, as well as for 
placing model-independent constraints on dark energy and modified gravity model building. 
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1 Introduction and Motivation 



As the accuracy of cosmological measurements continues to improve at a rapid pace, it is crucial to 
establish a robust and economic way to connect data with fundamental theory. In particle physics 
and condensed matter systems the framework of Effective Field Theory (EFT) has proven to be 
very successful at this endeavor (for reviews see [HE]). Recently, it has been argued that these 
methods can be useful in cosmology as applied to both inflation [31 IH [5] and late-time acceleration 
P El E] (see PH QUI EE] for related approaches). 

In [3], a covariant effective field theory for single field inflationary models was constructed, 
applicable both to the inflationary background and to perturbations about that background. This 
approach was generalized to the case of late time acceleration in [8] . A subtlety that arose in 
this generalization was the appearance of additional couplings and operators that involve matter 
and that give rise to modifications of gravity that can play an important role - appearing as a 
scalar-tensor theory in the regime of observational interest. Some observational constraints on the 
parameters of this type of EFT were derived in [12] , who showed through an attractor analysis that 
a large class of possible models could be eliminated by the data. 

This covariant, background- independent EFT approach has the goal of being able to describe 
all possible background cosmologies as well as perturbations. Because of this ambitious goal, it 
encompasses a smaller set of theories than an alternative EFT approach [HOE! [13] which fixes the 
background and attempts only to describe the theory of the perturbations. In particular, theories 
where higher derivative operators give order unity corrections to the background cosmology fall 
outside its domain of validity [$] . For inflation this includes models such as DBI inflation [H] , and 
for dark energy model it includes models like K-essence [15] and ghost condensates |16j . 

The alternative approach that has been advocated in the context of single field [I] and multi- 
field [5j inflation takes the background cosmology as given a priori and constructs the EFT of 
perturbations around that background. For inflation, much like the standard model vacuum, this 
makes sense as the inflationary background is not directly observable. It was noted in [3] that 
an inflationary background spontaneously breaks the time translation invariance of the underlying 
theory. This often makes it possible to capture the observational predictions of a class of inflationary 
models through the study of a single scalar mode - the Goldstone boson. The advantage of this 
approach is not only the simplicity of a single degree of freedom, but the clear explanation for the 
existence of universal relations between many of the parameters of the low-energy theory - such as 
those determining the sound speed and non-Gaussianity - that arise because the Goldstone must 
non-linearly realize the time translation invariance. 

In [6] , the authors applied this approach to the study of quintessence models to address issues 
of stability when models approach w < — 1. In this paper, we extend the work of [6] to allow more 
general couplings to matter, encompassing both dark energy and modified gravity models. We do 
so by generalizing the derivation of the theory of [I] to allow for the presence of an arbitrary matter 
sector that is assumed to satisfy the weak equivalence principle. We identify operators in the EFT 
that are interesting to study, and construct from those a general theory of scalar perturbations for 
dark energy and modified gravity. We also clarify some subtleties that arise in the formalism when 
matter fields are present; in particular, the scalar mode associated with the dark energy field is no 
longer the Goldstone boson of the broken time translation symmetry. 

With our theory in hand, we set about a program of surveying the interesting physical effects 
associated with different operators in the EFT. We calculate the contributions of each operator 
to a set of observables consisting of the speed of sound, the anisotropic shear stress, the effective 
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Newton's constant, Caldwell's w parameter that parameterizes modifications to gravity [T7], and 
other phenomenological parameterizations of modified gravity used in the literature. We conclude 
by summarizing the observational effects of different operators, and present a list of which operators 
are included in different dark energy models. 

Our work can be primarily viewed from two directions. From the effective field theory perspec- 
tive, it provides a description of all operators that can be present, and a unified framework for 
investigating their effects. From the perspective of comparing dark energy models to observational 
data, it provides a unified framework for describing the perturbative behavior of different mod- 
els. Each model can be matched to a certain set of operators, and with the framework developed 
here, the theory of linearized scalar perturbations within the context of the given model will be 
obtained. We anticipate that this formalism will be of use in constraining the size of different co- 
efficients through numerical simulations, both in model-dependent and independent manners. Our 
presentation of observational quantities is constructed with this in mind, and we have included full 
details of the equations of motion in both synchronous and Newtonian gauges in order to facilitate 
numerical analyses. 

While this paper was in preparation, Gubitosi et al. |13| published a very similar analysis. Their 
paper introduces some of the ideas of this paper, and as such, we have taken this into consideration 
when presenting our results. We refer readers to their paper for further details that have been 
omitted from this paper. Some aspects of our work that do not overlap with Gubitosi et al. are 
our explanation of the modified role of the Goldstone mode, the extent to which is valid to use a 
dark energy scalar field as a "clock" when it is a subdominant component of the Universe at high 
redshifts, the consideration of different operators in the EFT, the inclusion of spatial curvature 
throughout, and the casting of final results in synchronous gauge. Where there is overlap, our 
analyses are in agreement. 

The structure of the paper is as follows. In Section [2j we present an overview of the effective 
field theory construction, starting from the unitary gauge. From the resulting action, we identify 
how the background evolution is specified, and describe operators that affect only perturbations. 
We transform out of unitary gauge in Section [3j introducing the scalar field that describes dark 
energy perturbations. We briefly relate the resulting general description to a vanilla quintessence 
model, before investigating the meaning of the scalar perturbation field. In Section [4j we undertake 
a program of extracting information relevant to observations from our general theory, presenting 
the calculation of a variety of quantities. We also link our results to phenomenological functions and 
dark energy models. The paper concludes with a summary, and directions for future investigations. 

A number of detailed appendices accompany this paper. Appendix [A] details the application of 
the Stiickelberg technique for restoring full diffeomorphism invariance starting from a unitary gauge 
description. Appendix [B] gives some details of gauge transformation computations. We describe 
our conventions in momentum space in Appendix |C| and conclude by providing equations of motion 
in Appendix [Dj 

2 Summary of the Effective Field Theory Framework for Cosmo- 
logical Perturbations 

Our goal is to construct an EFT describing the perturbations about a cosmological background 
solution of a theory containing gravity, a single scalar field, and a matter sector which is assumed 
to obey the weak equivalence principle. The theory should include all possible couplings of the 
scalar to itself and to gravity, and will encompass both modified gravity and dark energy models. 
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Our approach is to follow the ideas of the EFT of Inflation [1] and the EFT of quintessence [6], 
generalizing to allow for arbitrary universal couplings to matter. In this section we describe the 
construction of the theory in unitary gauge. 

In constructing the theory, we assume that the matter sector obeys the Weak Equivalence 
Principle (WEP). Guided by this principle, matter fields will only appear in the matter action, and 
all matter fields will couple to the same metric. For our purposes, purposes, it is simplest to always 
work in the equivalent of the Jordan frame, where the metric variable g^ u is chosen to be the metric 
to which matter is minimally coupled. This choice is also convenient for comparing to observations, 
as the physical metric that we use for measurements must be the Jordan frame metricQ With this 
choice the matter action is an arbitrary functional S m \gui>,ip m ] °f the- Jordan frame metric and the 
matter fields ip m . 

We specialize to a choice of time coordinate for which the perturbation to the scalar field 
vanishes, that is, the class of unitary gauges in which the scalar field perturbation is "eaten" by 
the metric. The subgroup of the diffeomorphism symmetry group left unbroken by this gauge 
specialization consists of time-dependent spatial diffeomorphisms of the form t = t, x l = ^{t^x 3 ). 
We now use the tools of effective field theory to construct the most general action functional of 
a metric that is invariant under this group of residual symmetries. While the unitary gauge is 
not a useful gauge for investigating the behavior of perturbations, the allowed objects are easily 
recognized, and the action is remarkably straightforward in this gauge. We will transform later to 
more a more general class of gauges that makes the scalar perturbation mode manifest. 

The general unitary gauge action is given by 



Here, the quantities <5<? 00 , 5R, 5R^ U , and bK^y are perturbations (defined below) to the time-time 
component of the metric, the Ricci scalar, the Ricci tensor, and the extrinsic curvature tensor, 
respectively. The quantity C^ a \ is the Weyl tensor, and n M is the normal to surfaces of constant 
time. Each term has a coefficient that can be time-dependent, because the background solution 
breaks time translation symmetry. The final line contains S m , the matter action, which we leave 
arbitrary. The ellipsis indicates further terms that we have not included, at both quadratic and 
higher order in perturbations^] 

The perturbations of a quantity are defined by subtracting off the quantity's value on the 
background written in terms of the metric, normal, and functions of time. For example, we define 



1 This choice of frame removes some of the field redefinition freedom that was explored in Ref. [5]. 

2 The coefficients are largely chosen to coincide with the literature. However, we do caution the reader that some 
of our coefficients are different from Refs. [4] and [13]. In particular, our A(t) is defined differently (note that c(t) is 
the coefficient of Sg 00 instead of <? 00 ). There also exist differing choices of sign convention for the normal and K%j. 
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+ Mt){5R) 2 + X 2 (t)5R llu 5R^ + n(t)C^ aX C^ x + i2(t)e^ x C^P e C aXp g + ... 
+ m\(t)n"n v d tl 9™d v g m + m!(t)(<r + + . . .1 + S m [g^]. (2.1) 



SK^ = - K$ = + 3H(g flL , + n M rv). 



(2.2) 
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The ability to write the background components of operators in terms of dynamical variables is 
useful for maintaining symmetries in the action, and is unique to FRW backgrounds. 

The differences between the EFT of Inflation [4j action and (2.1) are threefold: 



• Our action contains a matter action S m . 

• Because we choose to describe the matter action in the Jordan frame, we now must include 
a function of time Q(t) multiplying the Einstein-Hilbert term. If we were to translate to an 
Einstein frame description, this function would parameterize the coupling of the scalar field 
to matter, as in a scalar tensor theory. 

• Terms such as 5R 2 must now be included in the action. In the EFT of inflation context, such 
terms could be eliminated using field redefinitions, but those field redefinitions are no longer 
allowed since they would take us out of the Jordan frame. 



Note that all of these differences stem from the inclusion of matter. 

The operators in the action split into two groups: those that contribute to the evolution of both 
the background and the perturbations, and those that only affect the perturbations. We will refer 
to the first group as the "background" terms (even though they do influence the perturbations). 
The background terms are given by the first line of the action (2.1), as well as the matter action. 



2.1 Background Evolution and Connection to ACDM 



We begin by deriving how the background terms in the action (|2.1|), given by 



So 



-9 



77) r 



n(t)R + A(t)-c(t)5g 



00 



+ S m [g, 



(2.3) 



relate to the background FRW solution, which is constrained by observations to be in close agree- 
ment with ACDM. Because the matter action couples universally to the metric and no other 
interactions are present, the matter stress-energy tensor T^( m ) is conserved: 



0. 



(2.4) 



Here, the (m) indicates a quantity associated with the matter fields. At the level of the background, 
we can treat the matter as a perfect fluid described by an average energy density p m and pressure 
P m , and so we have the usual continuity equation 



+ 3H[p m + 



0. 



(2.5) 



Here, and throughout the rest of the paper, overdots signify derivatives with respect to physical 
time. We do not use conformal time in this paper. 

A complete description of the background evolution entails knowledge of six functions of time: 



the scale factor a(t), the functions £l(t), A(t), and c(t) that appear in the action (2.3), and the 



matter variables p m (t), and P m (t). As we are constructing a theory of the perturbations, we require 
that the background evolution be specified a priori, and so all of these functions need to be chosen. 
We can use the Friedmann equations to eliminate two of the six functions, for which we will select 
A(t) and c(t). Typically, the functions a(t), p m (t) and P m (t) will be chosen to align closely with 
a ACDM evolution, but this does leave the function Q(t) free [181 E|- This extra function of time 
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describes, in the Jordan frame, a non-minimal coupling between the scalar field and the metric, 
and in the Einstein frame, a coupling of matter to the scalar field. 

The Einstein equations for the background evolution are obtained by varying the action (2.3) 
with respect to the metric. Using an FRW metric with spatial curvature ko, the Friedmann equa- 
tions are 



3rnkVL 



or 



n 



3H z + 2H + ^ + - + 2H- 



n 



Pr, 



A + 2c, 



-A - P m . 



For a given a(t), p m (t), P m (t) and £l(t), we can obtain c(t) and A(i) from 



(2.6) 
(2.7) 



c(t) 
A(t) 



Pm + P, 



m 2n 

m il 



■ k in HQ 

~^ + 20~ in 



Pn 



9 • kn £1 fi 
3H 2 + 2H+^ + - + 2H 
a z 



n n 

We may also write the Friedmann equations in a more conventional form as 



H + 



ko 



PQ+ Pn 



2mln 



H 



ko 

n2 



[pQ +P Q + P m + p Tl 



(2.8) 
(2.9) 

(2.10) 



where pg(t) and Pq{€) are the dark energy density and pressure respectively. The functions c(t) 
and A(t) can be related to these by 



(A _ PQ + P Q _ mln mlHh 
C[) ~ 2 2 + 2 ' 



A(t) = Pq- m 2 Q - 2mlH(l . 



(2.11) 



These equations describe how the background dark energy density and pressure evolve over the 
cosmological history. The terms proportional to and f2 are simply subtracting off the overcounting 
that comes from including the effective energy density and pressure from the non-minimal coupling 
U(t) in p Q and Pq. 

Combining the two Friedmann equations yields the background continuity equation. 



PQ + Pm = -3-ff [pq + Pq + P m + Pm] + SmlCl 



H 2 + % 



(2.12) 



Subtracting off the background continuity equation for matter (2.5) from this yields the continuity 
equation for dark energy. 



PQ 



-3H [p Q + P Q ] + 3mgO 



H 2 + 



ko 



(2.13) 



We now have an understanding of how the background evolution fixes the coefficients of the 
tadpole terms in the action. Before continuing to describe the remainder of the action, we pause 
to briefly discuss the importance of the background evolution. In previous work, a covariant EFT 
construction of the background evolution has been considered and perturbations of dark energy 
and matter have been examined 0IS]- It is well known that observations - such as standard candle 
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measurements to determine the expansion rate - taken alone cannot constrain the background 
evolution, as a high level of degeneracy in the theoretical parameters remains. The hope from 
the formulations [71 [8] is thus to attempt to differentiate regions of parameter space that are 
degenerate at the level of background evolution by using a perturbative analysis. However, by 
using a combination of theoretical and observational constraints some reduction in the allowed 
space of models is possible, as recently exemplified in [12] . where data in combination with an 
attractor analysis lead to some restrictions on models. Of note is that this analysis was unable to 
find a better fit to the observational data than the ACDM model. Thus, with the possible exception 
of early dark energy models that have a significantly modified cosmological evolution, using ACDM 
to approximately determine the behavior of the scale factor and matter energy density/pressure 
will provide a good description of the background. 

This helps motivate the EFT approach to the perturbations that we are taking here, which by 
construction is agnostic to the origin of the background evolution. In the absence of a separate 
analysis of how the background evolution of a model depends on its parameters, we cannot link 
the perturbative behavior of the model to its background evolution. In a sense, we give up on 
the background evolution and focus on the perturbations. We will evaluate later the observational 
significance of scalar perturbations in dark energy models. 



2.2 Effective Theory of Perturbations in the Dark Sector 

Having described the background terms, we now briefly describe the remaining terms in the general 



action (2.1). Each of these terms consists of a time-dependent coefficient and products of objects 
that vanish on the FRW background. 

The objects that we have available to us include 5g 00 , 5Kp U , 5Rp Up(7 (or Cp Upa ), SRp U , and 5R, 
where the perturbed quantities are defined above. Each operator is constructed from these objects 
pieced together with functions of time, g^ v ', n„, e fJ- " crX and derivatives Vm. By construction, the 
combination n^5Kp U is vanishing, and furthermore, the expression Rp U n^n u can be rewritten using 

R^nW = K 2 - K»K V » + V,(/V,/) - V v {n v V »n»). (2.14) 

Using these expressions and the ability to integrate by parts reduces the number of independent 
operators we need to consider. 

A particular operator can be classified by its order in the perturbed objects. As we are inter- 
ested only in linear perturbation theory, we will consider only operators that are quadratic in the 
perturbations. Even at quadratic order there are an infinite number of possible operators, because 
of the possibility of an arbitrary number of derivatives acting on the perturbed quantities. How- 
ever, each derivative increases the scaling dimension of the operator, and so operators with more 
derivatives become less and less relevant. 

It is convenient to split derivatives into space and time components using projection tensors. 
The projected space and time derivative operators are (5% + n^n u )Vp and n^n u Vp respectively. 
Operators involving derivatives (particularly time-projected derivatives) will often lead to the ap- 
pearance of higher order time derivatives in the equations of motion, which indicate the presence 
of new fields (which are often ghosts from the Ostrogradski instability). However, such terms are 
not problematic when treated within the context of an effective field theory [3] . 

Three-dimensional quantities such as the Riemann tensor of the spatial metric and contractions 
thereof do not need to be included in the EFT description, as such objects can be related to their 
four-dimensional counterparts through the Gauss-Codazzi relations [I], which involve objects that 
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have all been included in the EFT construction already. However, this does not preclude us from 
using the three-dimensional objects in our formalism, so long as we realize that they may not be 
independent from other operators. While from an EFT perspective, such terms are unnecessary, 
they can be useful for the purpose of matching to explicit theories of dark energy, and can also be 
useful for categorizing operators by separating out time derivatives from spatial derivatives in the 
four-dimensional gravitational invariants. 



A number of examples of possible operators are given in Eq. (2.1). Given so many operators, 
each of which associated with its own time-dependent coefficient, it is of interest to determine 
which operators should be retained and which can be neglected. For this question, we are informed 
by three perspectives. First, what operators are required to describe linear perturbation theory 
for existing dark energy models using this formalism? Second, what operators give rise to unique 
signatures of interesting physics? Third, how do the observational effects of an operator scale with 
energy? 

The third of these perspectives is of the most interest from the EFT point of view. Unfortu- 
nately, a complete analysis of this question is difficult to perform, as it is strongly dependent upon 
which operators dominate the kinetic dispersion relation of the scalar perturbations, and is also 
complicated by kinetic mixing with gravitational scalar perturbations. 

From our incomplete exploration of these questions, we propose that the following operators 
are most worth investigating (beyond the background operators): 

(5g 00 ) 2 , Sg 00 6K^, (SK^f, SK^SK^, 5g 00 5R^. (2.15) 

Here, R^ is the Ricci scalar of the spatial metric jjy, and 5R^ = R$' — 6/co- These operators 
are sufficient to describe linear perturbations about FRW cosmologies of very general single field 
dark energy models, including Horndeski's theory [19] I20j. In terms of physical signatures, they 
are sufficient to modify all cosmological signatures of GR that we have calculated. Finally, note 
that these operators are the leading order operators given their mass dimension, assuming a naive 
relativistic dispersion relation and ignoring kinetic mixing effects. The only operators at this 
order not included are (g^ + n ll n u )d ll g 00 d u g (>0 and n^n^d^g^d^g 00 , which contribute higher order 
derivatives to the equations of motion^} 

A number of these operators have been explored in the inflationary formalism, particularly with 
regards to correlation functions and stability. In particular, (5g 00 ) 2 , {5K^^) 2 and 5K^ U 5K V ' were 
investigated in Ref. [3], while (Sg 00 ) 2 and 5g 00 5K^^ were analyzed in Ref. [6]. The operators (5g 00 ) 2 
and 5g m 5K^^ were investigated in Ref. [13] . The operator 5g 00 5R^ has not been investigated 
before. 



3 Beyond Unitary Gauge 

In the previous section we showed that, in the unitary gauge, the EFT is specified by choosing a 
background solution (four free functions of time) and by choosing the time-dependent coefficients 
of a set of additional operators which only affect perturbations. We would like now to examine the 
predictions of this theory. To this end, is convenient to use a different representation of the theory 
in which the diffeomorhpism symmetry is restored and the dark energy field is explicit. 

3 Gubitosi et al. |13j suggested that the operator (g^ v + n ll n' y )d^g 00 d 1 yg 00 is of interest for describing Hofava- 
Lifshitz gravity. However, the equations of motion from this operator involve third-order time derivatives. While 
this operator may indeed be of interest, the analysis of operators involving higher time derivatives in the context a 
general theory is beyond the scope of this paper. 
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In this section we restore the full diffeomorphism symmetry in the action following Refs. [HE], 
by making use of the Stiickelberg trick, which introduces a scalar field tt, the dark energy field. We 
express the action to second order in perturbations, and briefly describe how this action can describe 
perturbations to vanilla quintessence models as an example. We then discuss some subtleties 
involved in the EFT construction in the present context where matter is present, and show that tt 
is not the Goldstone boson of broken time translations in this context. 



3.1 Restoring the Symmetry 

Given the action in unitary gauge, we can restore time diffeomorphism invariance through the use 
of the Stiickelberg technique, introducing a scalar field which nonlinearly realizes the symmetry. 
This procedure follows that given in the EFT of Inflation [3] , and we refer the reader to that paper 
for more details. 

The basic premise is to perform the following coordinate transformation on the action, 

t = t + Z°(x>*), x l = x\ (3.1) 

where we adopt the convention that Greek letters represent four-dimensional indices, while Ro- 
man indices represent three-dimensional spatial indices, and the zero index is reserved for the time 
component. The spatial coordinates can be left unchanged, as the action still possesses spatial dif- 
feomorphism invariance. Next, in the resulting action, replace £° everywhere by —tt(x ij '), effectively 
promoting £° to a dynamical field. The field tt is assumed to transform under time diffeomorphisms 
via tt — > tt — t;°(x^), so that t + tt is invariant. By construction, this will absorb any changes to 
the action caused by the change of coordinates. In effect, the action is now fully diffeomorphism 
invariant. The coordinates can be relabeled to remove the tildes. 

Using this procedure, explicit functions of time are modified according to 

f(t)^ f(t + 7r(xn) , (3.2) 

and will typically be Taylor-expanded in it. Operators that are not fully diffeomorphism invariant 
in four dimensions need to be transformed using the tensor transformation law. We describe this 
proceedure in detail in Appendix [Xj After this procedure has been performed, the final action is a 
functional of a metric and a scalar tt, and is invariant under time diffeomorphisms once again. 



Upon introducing the tt field to the action (2.1), including the operators in (2.15), we obtain 
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the following action. 

,2 



d x*J—g 



m. 



^n(t + ir)R + A(i + tt) 



c(t + vr) ( <5 5 uu - 2^ + 2^J 5 UU + 2V;7rcT - tT + 2_ ViTrVj-vr 



+ 



M|(t + 7T 



2 

Mf(i + 7T 



M|(t + 7T 



+ 



M|(t + 7T 



M 2 (t + 7T 



(«J B1 -M + ...) ! 



V 2 7T 

<5A^ + 3Hir +^r + 



V 2 tt 

+ 3HtT + — ^ + 



wrv + + ^VfeVj-Tr + ...)( ^ + + ^2 v « v ^ + • ■ ■ 



(5K\ + ...) 2 + 2[5K i t - H^Vjir + ...) (5K\ + HV^ + 



(^ 00 - 27T + ...)^ 3 ) + 



+ S m [g, 



1 1" i 



(3.3) 



The ellipsis in brackets reflect where we have truncated the expansions to quadratic order in per- 
turbations. The final ellipsis represents other quadratic terms that we have neglected, as well as 
terms beginning at third order in perturbations. Tildes are used to denote quantities associated 
with the spatial metric. 



3.2 Metrics and Gauges 



Although Eq. (3.3) is written in terms of an arbitrary metric perturbation, it becomes convenient 
to choose a gauge associated with metric perturbations. Given the present formalism where time 
and space have been treated on a separate footing, we find it convenient to work in synchronous 
gauge, where 5g 00 = g° l = 0. Well-known issues with this gauge will not arise for the purposes 
of our calculations. Furthermore, we envisage our equations being used in modified versions of 
software such as CAMB |21j . which uses synchronous gauge. 

We use the following metric to describe perturbations in synchronous gauge 

ds 2 = -dt 2 + a 2 {gij + h ij )dx i dx j , (3.4) 

where c/ij is the time-independent background spatial metric, and we let the spatial coordinate 
system be arbitrary. The perturbation hij is decomposed as 

hj = ~gij + [ViVj ~ |^V 2 ) r, , (3.5) 

where tildes indicate quantities associated with the spatial metric. The metric determinant y/—g 
is 

^g = a^g(l + f\. (3.6) 



11 



Following Ma and Bertschinger [22] , we use an alternative definition of r/ in momentum space, 
defined by 

r)(k,t) = -~[h(k,t) + 6fj(k,t)]. (3.7) 

We detail our conventions in momentum space in Appendix [C| 

It will be convenient to discuss some phenomena in Newtonian gauge, which avoids a number 
of issues associated with synchronous gauge. In Newtonian gauge, we use the following metric. 

ds 2 = -(1 + 2ij))dt 2 + a 2 (l - 2<j))g ij dx i dx j (3.8) 

We present the equations of motion for our theory in both gauges. For more information about our 
gauge choices, see Appendix [Bj 



3.3 The action for the dark energy field 



We will later need the general-gauge action (3.3) in order to calculate the stress-energy tensor. But 
for the moment, we are interested in the dynamics of the ir field. Expanding (3.3) to second order 
in perturbations in synchronous gauge, the it action becomes 



,4 % i ~ 

d xa \J g 



c TT 



,2 9°- 



VivrV^TT + 3Hc 



m. 



+ 2M|tt 2 + Mf-k 



+ 3Hir + 



4 

vV 



nip 



2 



h ■ V 2 7T 

- + 3H7T+ 



h'; 



V 1 



+ Ml(t)H 2 ^- VivrVj-vr - M 2 n5R {3) + 



where we have used the background continuity equation (2.13) in the form 

,2 



^OR (0) + A-2c-6Hc = 0. 



(3.9) 



(3.10) 



to cancel tadpole terms and simplify the action. Here, represents the background Ricci scalar 
as a function of time, and SR = R — R^ , which is given by 



fen 



1 



2 ~ 



2 ~ 



SR = AHh + h- + ( 2k V z rj - -V z h + -V 2 V 2 r/ 



(3.11) 



in synchronous gauge. Note that there is a kinetic mixing between the gravitational perturbations 
and 7r, arising from a number of different operators. Because of this, we cannot read the speed of 
sound directly off this action. We present the equations of motion from this action in Appendix |D} 



3.4 Example: The EFT of Quintessence 

In order to elucidate the ir field formalism, here we write a simple quintessence model using our 
language. This example follows the work of Creminelli et al |6j. 
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Consider the following minimally coupled quintessence action. 



nil 



-R 



{Vcj>f 



U(cj>) 



(3.12) 



Assume some cosmological background solution <f> = 4>o(t) on an FRW metric, and ignore metric 
perturbations. The solution <pQ obeys the background equation of motion 



□0o - U'{<fo) = 0. 

The background pressure and energy density are given by 



(3.13) 



(3.14) 



Now, expand in perturbations as (f> = (j)o(t + ir) = (j>o(t) + 4>qtt + ^4>qK 2 + • • •• I n terms of the usual 
quintessence perturbation <p = 00 + Scj), we see that 5(f) = 0O 71 " to leading order in perturbations. 
Using this expansion in the action (3.12), we obtain 

d A x 



-9 



m P D ^„00i2 



2 J R-- 5 uu ^ + (n0 o ) 



ivr + 



12 2 
b TT 



;[V(0o^] S 



u((t> ) - u'(M 



1 ■ 



2^2 
0^ 



(3.15) 



after performing an integration by parts. The background equation of motion can be used to 
eliminate the terms in the large round brackets. After expanding derivatives, integrating by parts, 
using the time derivative of Eq. (3.13) to simplify terms, writing g 00 = —1 + Sg 00 and using the 
background expressions for energy density and pressure, the resulting action is 



S,h 



d A XyJ—g 



m. 



-R + P6 



P4> + P 4 



<t> s: 00 



Sg 



(Vvr) 2 + 



(3.16) 



By comparing the background terms to Eq. (2.3), we conclude that 
A(t) = P q 



U ((/>), c(t) 



2 ' 



i 



(3.17) 



2 ' 2 

for this model. The remaining terms then agree with our general expansion, Eq. (3.9), with all 
other operators identically vanishing. 



3.5 The EFT Construction in the Presence of Matter: Some Theoretical Con- 
siderations 

As we described above, the basic ideas of the original construction of the EFT of Inflation [3] are 
as follows. First, specialize to the unitary gauge associated with the inflaton field, so that the 
perturbation to the inflaton field vanishes, and the only remaining field is the metric. Second, 
write down the most general action that respects the subgroup of the diffeomorphism symmetry 
group consistent with that gauge condition. Third, use the Stiickelberg technique of introducing a 
new dynamical field, ir, to restore full diffeomorphism invariance. The new dynamical field turns 
out to be the Goldstone boson associated with the breaking of time translation invariance by the 
background solution. 

Consider now a more complicated context such as multifield inflation [5] , or the case considered 
in this paper, where matter and a dark energy field are both present. In these cases, a number of 
questions about the EFT construction arise: 
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• Which unitary gauge do we use? For example, if we have two scalar fields (ft and ip, do we 
use the gauge for which 5(ft = 0, or the gauge where 5ift = 0? If one field is the dominant 
component of the energy density in the Universe, must we use that field to define the gauge? 

• When we use the Stiickelberg technique to restore full diffeomorphism invariance, is the 
scalar field n that is used still the Goldstone boson associated with broken time translation 
invariance? 

The answers to these questions are, first, that at a classical level it is legitimate to use unitary 
gauge associated with any field, even a subdominant one, to define the theory. The only compli- 
cation is that the usual EFT rules for estimating the coefficients of operators in terms of a cutoff 
scale need not apply if one uses a subdominant field, as we discuss further below. Second, the field 
7r that arises from the Stiickelberg technique is no longer the Goldstone mode in general. 

To explain these issues, we will specialize temporarily to a simple model system of a particle 
moving on a iV-dimensional manifold whose action is 

S[x A (r)} = J ^{^[x C (r)]^^-y[x c (r)]}, (3.18) 

where x A = (x 1 , . . . ,x N ). This system is not invariant under time reparameterizations. However 
we can rewrite it in a way that is by introducing the one dimensional metric dr 2 = e 2a ^dt 2 , in 
terms of which the action becomes 

S[x A (t),a(t)} = J dt\^e-^g AB [x c {t)]^^-e^V[x c {t)\Y (3.19) 
This form of the action is explicitly invariant under the time reparameterizations 

t->i=i(t), x A (t) -»• x A (t) = x A (t(t)), a{t) -»• a(t) = a(t(t)) + In (J^J . (3.20) 

Consider first the status of the Goldstone mode in this model system. Specialize for simplicity 
to the gauge a = 0. Suppose that we have a background solution 

x A (t)=x^ A (t), (3.21) 

and we are interested in the theory of the perturbations 5x A (t). The Goldstone mode ir associated 
with the broken time translation symmetry is a coordinate on field space defined by 

5x A {t) = x (0)A {t + tt(*)) - x (0) A {t) , (3.22) 

for 1 < A < N. To get a complete set of coordinates on field space one has to add N — 1 other 
coordinates, say, £ 2 , . . . ,£ N . The Goldstone mode has the property that it = constant, £ A = for 
2 < A < N, is always a solution. The Goldstone mode is only defined when we are given a complete 
background solution. 

Consider next the Stiickelberg technique. This trick can be defined in very general contexts, 
when one does not specify a background solution, and so it cannot involve the Goldstone mode in 
general. As an example, let us fix the gauge in our model by demanding that 

x 1 (t) = f(t), (3.23) 
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where f(t) is some fixed function that we have picked (it could be a background solution for x 1 ). 



The action (3.19) becomes 



S[x r (t),a(t)]= J ^|I e -VE[/(t),x e (t)]^^ + e-Vi[/W,x e (t)]^| 

+ J dt^e-«g u [f(t),x e (t)] (|) 2 -e^)F[/(t) )a; (t)]|. (3.24) 

Here Greek letters T run over 2, . . . , N, so x A = (x 1 , x r ). Note that now the field a is no longer 
pure gauge, it is physical. 

We now want to restore time reparameterization invariance. Following the procedure outlined 
in [3], we make a coordinate transformation t — > t = t + £(t), rewrite all tensor components in terms 
of their transformed versions, replace £ everywhere by w, and then drop the tildes. This gives the 
action 

*M>Mt)Mt)] - J * + S«» + 

+ y<it|V"(') 911 (^/(t + fW)) 2 -e' ,(,, vl , (3.25) 
where the dependence of the metric and potential on [f(t + vr(t)), x @ (t)] has been suppressed. By 



inspection we see that this is the same as the original action (3.19), if we identify fit + vr(t)) with 
x l {t). Thus the Stiickelberg trick does what we expect in this context: at a classical level, any 
unitary gauge can be used, and one can freely transform back and forth between the unitary gauge 
and general gauge representations of the theory. 

Now suppose that one is given a background solution x^ A (t). Then, for each A, one can define 
a unitary gauge by x A (t) = x^°> A (t). Returning to a general gauge via the Stiickelberg trick, for a 
given fixed A, will involve a field ir A defined by 

x A {t) = x {0)A [t + 7r A (t)]. (3.26) 

The quantities tt 1 , . . . ,ir are a set of coordinates on the manifold that can be used instead of 
x 1 , . . . , x N . In these coordinates the Goldstone mode is the diagonal excitation tt 1 = tt 2 . . . = tt n = 

TT. 

Finally, consider a cosmological model with one dominant field and one subdominant field. If 
one uses the unitary gauge associated with the dominant field (i.e., uses that field as a clock), 
then the subdominant field can be completely neglected, and the situation reduces to the single 
field EFT case. Here, the standard EFT rules for the scaling of coefficients of operators with the 
cutoff of the EFT will apply. On the other hand, suppose one uses the unitary gauge associated 
with the subdominant field. At a purely classical level, this is a legitimate thing to do. However, 
because one is using a gauge that is poorly adapted to the dominant physics, the usual EFT scaling 
rules for coefficients of operators will no longer apply. Specifically, there is a new dimensionless 
quantity (roughly the ratio of energy densities in the dominant and subdominant fields^]) on which 
the coefficients can depend. Therefore a subdominant field is not a very useful choice for a clock 
to define the gauge. 



4 A more careful analysis suggests that the ratio of sum P + p of each component may be the relevant quantity to 
consider. Note that this sum is responsible for breaking the time translation symmetry in the acceleration equation, 

lOoi. 
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In the application of this paper, the energy densities of dark matter and dark energy are 
comparable today, and so our construction will be valid at low redshifts, with the usual EFT 
scaling rules being applicable. However, at high redshifts, the dark energy field which we are using 
as a clock will be a subdominant component. This suggests that care must be taken when naively 
applying the EFT analysis to observations at high redshifts, since the breaking of time translation 
invariance at those redshifts is primarily due to dark matter and not dark energy. Understanding 
the EFT approach at high redshifts when dark energy is subdominant represents work in progress. 

4 Linking the Behavior of Perturbations to Observations 

We now consider the full system of perturbations around the FRW background. We use the 
linearizes equations of motion for scalar perturbations to calculate the speed of sound, Poisson 
equation, anisotropic stress, effective Newton's constant, and Caldwell's w parameter. We associate 
our constructions with phenomenological functions appearing in the literature, and conclude by 
highlighting the observational features of different models. 

Some of these quantities have been calculated by Gubitosi et al. [13J, using the operators M2 
and M\. Their calculations were mostly performed at the level of the action, by explicitly demixing 
the field content through field redefinitions. Our approach is related, but works entirely at the 
level of the equations of motion. The benefit to this is that we can compute physical quantities 
in terms of equations of motion with arbitrary coefficients. We then present the contribution 
that different operators make to these coefficients, allowing these quantities to be calculated for 
whichever combination of operators is desired. 

4.1 The Role of Perturbations 

As mentioned earlier, it is impossible to discriminate between different dark energy models based on 
their background evolution alone; we need to investigate how different models impact the evolution 
of perturbations, particularly the growth of large scale structure. Typically, this is influenced by 
the cosmological history, with little impact from dark energy perturbations. Once dark energy 
comes to dominate the universe (around z ~ 2), the perturbations may play a more prominent role, 
but are still unlikely to have a sufficiently noticeable impact on structure formation, with possible 
exceptions for small sound speeds. Unlike inflation, where the scalar perturbations set the initial 
conditions for future evolution, we do not expect to measure a correlation function associated with a 
dark energy scalar perturbation. It is thus worth asking why we are interested in the perturbations 
described by the ir field? 

The answer is twofold. Firstly, understanding the dynamics of the perturbations is critical 
to controlling the stability of the theory. Even when dark energy is subdominant, an instability 
on a sufficiently small timescale can have catastrophic consequences. The issue of stability for 
quintessence has been investigated in some detail in Ref. [6]. 

Secondly, and more interestingly, when there exists a kinetic mixing between the scalar field 
and the gravitational scalar perturbations, the gravitational perturbations have their propagation 
modified. As ir perturbations in the form of it source the gravitational scalar perturbations, the ir 
equation of motion, which contains second derivatives of the gravitational scalars, must be substi- 
tuted in. This can cause significant modification to the dynamics of the gravitational perturbations. 
As these perturbations are also sourced by and react upon the matter perturbations, the behavior 
and growth of matter over- and under-densities are modified in this situation. Even when dark 
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energy is a subdominant component of the universe, it can have a nontrivial effect on the growth 
of structure through this modification of the gravitational dynamics. 

Therefore, this formalism is most powerful for constraining theories involving modified gravita- 
tional dynamics. In particular, we envisage that the size of some coefficients should be able to be 
constrained based on the drastic departure from general relativity that they entail. Conversely, for 
less "exciting" theories such as minimally coupled quintessence//c-essence, an analysis of large scale 
structure formation will provide information about the background evolution of the universe more 
so than the size of various coefficients. This is not particularly surprising, given the small number 
of operators necessary to describe minimally coupled quintessence and /c-essence models. Never- 
theless, the ability to constrain A(t), c(t) and Q(t) is important, as these functions can describe 
basic deviations to the cosmological evolution as compared to ACDM. 

This discussion highlights why the present formalism is an important extension of our previous 
work 0|8], in which the regime of validity demanded that the operators corresponding to f2, c and 
A dominated the perturbative behavior, while other operators played a subdominant role. 



4.2 Speed of Sound 



To calculate the speed of sound, we need to take the determinant of the kinetic matrix that mixes 
7r, h and fj. Setting the determinant to zero yields the dispersion relationship for the fields. We 
can form this matrix from the equations of motion, given in Appendix [D] 

After replacing time derivatives with —icj, write the kinetic terms (neglecting dissipative terms) 
from the equation of motion for tt, the time-time component of the Einstein equation, and the 
space-space trace component of the Einstein equation as 



^ ln,n ln,h 77T,7j \ ( K ^ 

h 

\ V J 



7tt,n ltt,h 7tt,?? 
\ 'Jss.n "JsSjh 1ss,fj J 



(4.1) 



The first index is intended to indicate the equation of motion the component comes from. Note 
that we are using synchronous gauge for this calculation. 

The contributions to this kinetic matrix from various operators are detailed in Table [l| The 
contribution to the kinetic matrix from the background operators is worth paying attention to, as 



this fixes the normalization and overall sign for each row. In particular, we use Eqs. (D.l), (D.13) 



and (D.15) to extract the background contribution. The ir equation of motion has been left largely 
unchanged, while the gravitational equations of motion have been written as = —G^v + T% . 

As an example, the determinant of the kinetic matrix from the background operators alone 
yields 



rtlj 



C + 



UJ 



c + 








(4.2) 



with a meaningless multiplicative constant, indicating that there is only one propagating degree of 
freedom, with speed of sound unity. 

The dispersion relationship arising from the operators M2 and M3 is more complicated than 
a relativistic scaling, due to the presence of fourth order spatial derivatives in the action. The 
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Operator 



Kinetic Matrix Contribution 



Speed of Sound c 2 



Q(t),c(t),A(t) 



2 



Q 



7^ 



2 fc z 

3 a* 
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, 3 2 Si 2 
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6g m SK^ 
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Mi \ fc 2 
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'1 2?T 
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4 m 2 n 
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2 o 



r — 3H ^2 

3// 4 







V (3^ + 2f )g -a> a /2 J 
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SK^SK" 



Mi 



H_ k A 
2 a 2 " 
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-w 2 /6 y 



(See text) 



(See text) 



M 2 (t) 



Sg 00 5R i3) 



M 



' -2 (3# + 2|f) fc 2 \ 

-4fc 2 
V / 



c + |m 2 f - ga 2 (3H + a t )M 2 + a 2 M 2 ^ + 2c-° 



! 2 S2 



c+ | m?) «i + | a 2M 2 «i + 2 C ^ 



Table 1: The contribution to the kinetic matrix from various operators. The last column gives 
the speed of sound for operators combined with the background operators only. 



dispersion relation from M2 combined with background operators is 



a; = 



2 m 2 Q 2 m 2 ^ ^ 4 mgfi f7 



4'"'0 Q 2^7 



a 2 „ 1 3^2 SI 2 1 3cA? 2 2 

(4.3) 



a 4 r , 3_2£i 2 . , 3 egg 



while the dispersion relation from M3 combined with background operators is 

2 *« M|/2 + |§ 
y i_ 

~ a 4 _ , 3 rrj 2f2£ , 3cM| 

c + 4 m o n + 2^ 

,9 3 m 2S7 2 , la ,/V>2\f2 / H 2 , m )Tf 2 , 3 ir,(>2fi , 3 A >2 Q 2 , 3 cM| 3 (g+ff 2 )M 3 4 3 ™^ 

(4.4) 
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It is possible, in principle, to apply a unitary transformation to the fields ir, h and fj to demix 
the fields at the level of the action, as described in Ref. [13] . The kinetic term in the action would 
then look like 



inetic 



/ 



d 4 Xy 



-g 



urn 



i'2 



(4.5) 



where tt is the propagating field. This form is of interest, because it facilitates investigations of 
stability. The speed of sound can be read off from this action as c 2 = g(t)/ f(t). However, there is 
an ambiguity in relating f(t) and g{t) to the coefficients of u 2 and k 2 /a 2 in Eq. (4.2) and similar 
dispersion results, as / — > fh(t),g — > g/h(t) will yield the same dispersion relationship. We can 
fix this ambiguity by taking the limit where fi — >• 1 and only the background operators are turned 
on, so that we know that f(t) = g(t) = c(t) from the action (3.9). When the dispersion relation 
changes, the kinetic action for the propagating degree of freedom can be constructed in a similar 
way, although this will require the addition of new operators to the action (4.5). In Table [TJ the 
speed of sound is written as g(t)/ f(t). 

The speed of sound is intimately linked to questions of stability. Here, we comment only very 
briefly on some salient points. By requiring that the theory is ghost-free, the coefficient f(t) of 7r 2 



in Eq. (4.5) must be positive. For just the background operators, this requires 



,^2 



PQ+Pq 



+ 



m 2 Htl 



>0, 



(4.6) 



which indicates that the equation of state of dark energy can become phantom-like (wq < 1) without 
causing perturbations to become ghosts, for a suitable choice of Q(t). The presence of further 
operators creates more mechanisms to avoid ghosts when the equation of state becomes phantom- 
like. If the numerator g(t) becomes negative, then a gradient instability arises. Furthermore, a 
Jeans instability can also arise from certain operators. These issues were explored in detail in Ref. 

El. 



4.3 Poisson Equation and Anisotropic Shear Stress 

The Poisson equation combines the time-time and time-space components of the gravitational 
equations of motion in Newtonian gauge. Breaking the gravitational equation of motion into pieces 



following Eq. (D.6), the Y mode of the time-time component can be written as 

'2k 2 - 6k 



+ 6H(<P N + H^) 



s p n + r 



,0 (Q) 



Y 



while the Y{ mode of the time-space component is given by 



mln(t)2k 



l N + Hi; N 



(p m + P m )v N + T 



,0 (Q) 



y, 



Substituting the second of these equations into the first, we obtain the Poisson equation. 



2 2k 2 -6k N 
m U g <p 



-5 A + T° {Q) 



-0 (Q) 



Y k 



Yi 



Here, we have defined 



p m A = 5p N + 3 ^{p m + Pm)v N , 



(4.7) 



(4.8) 



(4.9) 



(4.10) 
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Table 2: The contributions to the Poisson equation and anisotropic shear stress from various 
operators. Note that for the Poisson equation contribution from the background operators, the 



background continuity equation (2.13) has been used to simplify the result. 



as the rest-frame density perturbation of matter. 

The anisotropic shear stress is given by the Y l - mode of the spatial components of the gravita- 
tional equation of motion. We express it in Newtonian gauge. 

k 2 



m 2 n- 2 - 



p m n + t 



i (Q) 



Y': 



(4.11) 



In order to detail the contributions to the Poisson equation and anisotropic shear stress from 
our operators, we define Ap for an operator by its contribution to the Poisson equation by 

2k 2 



6&0 



N 



Qiu a 2 



-p m A + . . . + Ap. 



(4.12) 



[cf> N -i; N )=P m U + ... + Ac 



(4.13) 



We similarly define Ag for the anisotropic shear stress as 

k 2 
a 1 

In Table [2j we present the contribution to the Poisson equation and anisotropic shear stress from 
all operators we considered. These quantities were extracted from the equations of motion, detailed 
in Appendix [D} 



4.4 Newtonian Limit 

In order to calculate the effective Newtonian constant, we need to take the Newtonian limit of our 
equations of motion. The Poisson equation, the anisotropic stress equation, and the equation of 
motion for tt contain the information we need. We follow the method described in Ref. [23] closely. 

In the Newtonian limit, we look at modes that are deep within the horizoij^J k 2 /a 2 3> H 2 . 
We also make the quasistatic approximation, discarding time derivatives of the fields. Finally, we 
assume that matter gives rise to no anisotropic shear stress. 



De Felice et al. [23] point out that if the speed of sound is much less than one, then the sound horizon should be 
used instead of the Hubble horizon, which may limit the regime of validity of this approximation substantially. 
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Table 3: The contributions to the Poisson equation, anisotropic shear stress, and tt equations of 
motion in the Newtonian limit. 



To include the effect of possible perturbative terms, we write the Poisson equation, anisotropic 
stress equation, and it equation of motion in this limit with arbitrary coefficients: 



A: 2 



+ A 



k 



r7T 



N 



,A. 



BnP N + B 2 cj> N + B 3 ir N ~ , 



Ci 



7J N ~ 



0. 



(4.14a) 
(4.14b) 

(4.14c) 



We keep the "mass" term M 2 in the 7r equation of motion; whether or not this term is suppressed 
will depend on how it compares to the mode in question. Note that M 2 actually has mass dimension 
6, because the tt field is not canonically normalized. In Table [3j we present the contribution to 
each coefficient from the operators we have considered. The overall sign and normalization for 
the coefficients are given as follows. The Ai coefficients can be read off from — Ap, while the Bi 
coefficients can be read off from — Aga 2 /^m^O, (the normalization is chosen for simplicity). For the 



background contribution, the LHS of Eqs. (4.9) and (4.11) need to also be included. Contributions 



to the Ci coefficients and M 2 can be read off from A,,-, the contribution of individual operators to 
the 7r equation of motion, which is presented in Appendix |Dj 

To obtain the effective Newton constant, we want to form the equation 

k 2 



-47rG cff/ o m A. 



(4.15) 



Take Eq. ( |4.14b[ ), solve it for 
(414c) for V v 



and substitute this into Eqs. (4.14a) and (4.14c). Then, solve 



substitute this into (4.14a), and solve for G e g. 

B 2 C 3 - dB 3 + ^B 2 M 2 



4vrG eff 



A 1 B 3 C 2 - A 2 B 2 C 2 + A 2 Ci5i - AiBiCs - ^A X B X M 2 



(4.16) 
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2m 2 1 fi ^ 40 c + 3mgn 2 + 4fiM 2 - 12VHM 2 + 2Q.M 2 g + 3ffM|0 / 




1 / rngfi 2 - + 2d t )M! + 2H Hdt{&2) \ 
11 
2m 2 fi I 4 rj c + 3m 2f22 + 4Q.^M 2 - 4M 2 H 2 Q + 2M| 0,^ - 4M 2 QH - Q(3H + 2S t )M| - 2H^d t (M 2 ) J 


2 


1 ( mid 2 4a 2 Q(3H + d t )M 2 4a 2 M 2n '" 8c 8°"Tm 2 4 \ 

I 1 —|— m m 

zm ii y 40e + 3mg0 2 + 4Q^M 2 — 4a 2 Q(SH + d t )M 2 + 4a 2 M 2 ^- + 8c ^" + 8 M ^2 J 


Operator 




General Relativity 





Q0),c(*),A(t) 


™2 n 2 
c + m 2 g + 2^M 2 




1 o 2 

c + m ?M + T2 M ' 2 


Sg 00 SK» 


2 si 2 ^3 n 




c + ™§ i + !^ 2 " 1 (H + dt)M? + Ml A 




m 2 g_|M 2 i/g 

c + ^ 2 i + f^^ 2 + ^ 2 2 ^ - 3A? 2 ^ + :\/^'//;; 




mgi - |j(2i/ + 9 t )M| + + a t )M| 
c + mg^ + f^M 2 - M 2 (H 2 + H) + M 2 ^ - ^(M 2 ) - gg(tf + a t )M| 


2 y 


m oS ~ §a 2 (6H + 2<9 t )M 2 


Table 4: 


Effective Newtonian constants and w parameters for models involving only the back- 



ground operators and individual further operators. The values for general relativity are included 
for comparison. The mass scale M 2 is not written out explicitly, but can be read from Table |3j 
Note that for most operators, the contribution to w is nonzero only for Q ^ 0. 
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We can also use Eqs. (4.14b) and (4.14b) to eliminate 7^ and solve for w = ip N /4> N — 1, where 



w is defined by Caldwell el al. 



|17j as the cosmological analog of Eddington's PPN 7 parameter 



w 



B2C3 — C1B3 + j^B 2 M 2 
B 3 C 2 - B 1 C 3 - %B X M 2 



(4.17) 



This result for w is calculated from a quadratic action, and so we stress that nonlinear effects such 
as the Vainshtein [24|. 125] [26] and chameleon |27l [28] mechanisms are not present. Thus, solar 
system constraints on 7 cannot be used to constrain w. However, cosmological observations may 
be used to do so; Caldwell et al. [T7] place a bound of w = 0.02 ± 0.07 (2a) on tens of kiloparsec 
scales, and w = 0.03 ± 0.10 (2a) on hundreds of kiloparsec scales. Unfortunately, no useful limits 
were obtained on Megaparsec scales. 

We present the effective Newtonian constants and w parameters for models involving just the 
background and individual further operators in Table [4] Just as with the speed of sound, more 
general results can be calculated by combining the effects of multiple operators. 



4.5 Phenomenological Functions 

In order to compare observational data with modification to general relativity, Bean and Tang- 
matitham |29] introduced two phenomenological functions Q and R which perturbed the Poisson 
equation and anisotropic shear stress. Here, we construct their functions from our formalism. Other 
phenomenological parameterizations have also been proposed; see |29] and references therein. 

As we will be performing comparisons to general relativity in this section, we need to decide 
what form Newton's constant G will take in the comparisons. For simplicity, we use rn^fl = 1/8ttG 
here. However, it may be preferable to absorb O into the phenomenological functions, and it will 
be important to decide which scales to compare to the Newton's constant measured today in order 
to perform a thorough analysis. 

In general relativity, the Poisson equation is 

(k 2 - 3k )(j) N = - 47rGo 2 yO m A , (4.18) 
where A is defined in Eq. (4.10), while the anisotropic shear stress is given b>|^] 



^_^ = 8vrGa 2 ^. (4.19) 

Bean and Tangmatitham introduce two phenomenological functions Q and R to describe deviations 
from general relativity by 

{k 2 - 3k )(j) N = - 47tGQa 2 p m A , (4.20) 

and 

<t? -R^ N = 8 ^Ga 2 ^. (4.21) 



6 An alternative representation of the matter anisotropic shear (such as used by Ma and Bertschinger [22]) is given 
by a = 2P m n/3(p m + 
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We may write our Poisson equation (4.9) in this form as 



c 2 - 3^)0^ 



1 + 



T 



(Q) 



3H o (Q) 



and thus identify 



Q = 1 



(4.22) 



(4.23) 



where the sum is intended to indicate the combined contribution over all operators. Similarly, we 



can write our equation for the anisotropic shear stress (4.11) as 



mlp, k 2 



leading to 



1 a 2 

R=1 + ^m 2 tW 2 ^ As - 



(4.24) 



(4.25) 



For a typical dark energy model, the function Q is expected to remain close to unity until the 
very late stages of the universe, as the energy density in dark energy perturbations will be dwarfed 
by that of the matter perturbations. However, the effect of a non-minimal coupling 0, will change 
the evolutionary dynamics from the time variation of the Newtonian "constant" . The R parameter 
is more interesting, and displays the possibility of deviating from unity much earlier, as it is the 
ratio of the dark energy perturbations to the gravitational perturbations which is important. This 
particular combination is probed by weak lensing surveys. 

Bean and Tangmatitham found that a constant Q and R were consistent with unity at the 
95% confidence level. However, for Q and R that became important in the late universe, much 
more freedom was available, particularly on large scales. At the time that work was performed, 
observational data was insufficient to lift a degeneracy between Q and R. It is hoped that newer 
data will be able to lift this degeneracy. 



4.6 Confronting Models with Observations 

Having calculated the contribution from various operators in our EFT action to different observ- 
ables, we can now combine this with our knowledge of different dark energy models. Table [5] 
summarizes the operators required to describe various models. 

We have calculated three main quantities. The first is the speed of sound, which is closely related 
to stability considerations, and also to the clumping of dark energy. The dispersion relationship 
may also be modified away from the relativistic oo = c s k. The second is the modified Poisson 
equation in the Newtonian limit, which is important in describing the growth of structure. We 
describe the modified Poisson equation in terms of the effective Newtonian constant. The third is 
the effective anisotropic shear stress, described by Caldwell's w parameter, which appears in weak 
lensing equations. Each operator gives a different contribution to each of these quantities. 

Before summarizing the effect of each operator on these quantities, we should begin by un- 
derstanding what is predicted by the ACDM and basic quintessence models, so we know what we 
are comparing to. In ACDM, dark energy exists in the form of a cosmological constant. There is 
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Operator 


n 


A 


c 




Mf 


Mf 




M 2 


Model 


R 




5<? 00 


(V ) 2 


5g 00 6K^ 








ACDM 


1 


/ 





- 


- 


- 


- 


- 


Quintessence 


1// 


/ 


/ 












F(R) 


/ 


/ 















fc-essence 


1 * 


/ 


/ 


/ 










Galileon 

Kinetic Braiding 


1// 


/ 


/ 


/ 


/ 








DGP 


/ 


/t 


/t 


/t 


/ 








Ghost Condensate 


1// 


/ 









/ 


/ 




Horndeski 


/ 


/ 


/ 


/ 


/ 


/t 


/t 





Table 5: A list of operators required to describe various different models. The correspondence to 
various models is given in Ref. [T3] . The correspondence to Horndeski's general theory [p5] is given 
inRef. [20]. 

■S Operator is necessary 

Operator is not included 
1, Coefficient is unity or exactly vanishing 

1// Minimally and non- minimally coupled versions of this model exist 

f Coefficients marked with a dagger are linearly related to other coeffi- 
cients in that model by numerical coefficients 

no speed of sound and no clumping of dark energy. The effective Newtonian constant is exactly 
the Newtonian constant Gn, and anisotropic shear stresses arise only from photons and neutrinos, 
which are of negligible importance in the late universe. 

In contrast, a minimally coupled quintessence model has nonzero A(i) and c(t), and a speed 
of sound equal to the speed of light. The effective Newtonian constant is once again Gn (in the 
Newtonian limit). There are no additional sources of anisotropic shear stress. Now, as additional 
operators are turned on, we begin to see new effects appearing. We summarize these effects in 
Table 1 

By combining the data in Tables [5] and [6j we can predict which models will manifest different 
effects, and use this information to discriminate between different classes of models. It is hoped that 
with the help of numerical simulations, the coefficients of different operators can be constrained, 
leading to an overall picture of which classes of models yield a viable description of the growth of 
large scale structure. 
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Operator 


fi(t) 


Ml 


Mf 




M| 


M 2 


Observable 


R 


(dg 00 ) 2 


Sg°°8K^ 


(«^) 2 


SK^K^ 


6g 00 5RW 


Speed of Sound 


1 


/ 


/ 


/+ fc 4 


/+ fc 4 


/t 


Effective Newtonian Constant 


/ 




/ 


/t + k 2 


/+ fc 2 


/ 


Caldwell's w parameter 


/ 




/t 


/t 


/ 





Table 6: The effect of various operators on physical quantities. 

/ Operator modifies this quantity 

/ f Operator modifies this quantity when 0(f) 7^ 

Operator does not modify this quantity 
1 Speed of sound is unity 
/+ k n Operator introduces a new scale dependence 



5 Conclusions and Outlook 

We have proposed an effective field theory construction to describe dark energy and modified gravity 



models in the late universe. The construction of the action (2.1) is performed in a unitary gauge, 
and consists of operators written in terms of deviations from an FRW background of quantities that 
are invariant under spatial coordinate transformations. Three of these EFT terms, along with the 
matter action, contribute to the background cosmological evolution. A number of further operators 
deemed interesting were identified. 

Explicit scalar perturbations in our model were introduced in Section [3j where we also raised 
theoretical issues associated with the use of this construction at high redshifts. A systematic 
investigation into the physical effects of various operators was undertaken in Section |4j where we 
highlighted the importance of perturbative effects in constraining dark energy and modified gravity 
models. 

On the theoretical side, further work towards understanding the details of the formalism at high 
redshifts is necessary. Following this, a detailed investigation into which operators are dominant 
at different energy scales and cosmological epochs would also be of interest. These calculations are 
complicated by the presence of both the matter fields and the non- minimal coupling f2, and thus 
results from inflationary work using this formalism are non-trivial to translate. 

On the numerical side, this formalism represents an exciting opportunity to place model- 
independent constraints on dark energy and modified gravity models. Each individual operator 
leads to a characteristic physical signature, and follow-up work to constrain the behavior of the 
coefficients of the operators investigated here using observational data is underway [30] . 

This formalism will also be of use in investigating explicit models of dark energy, by linking the 
behavior of perturbations in specific models to the parameters which also control their background 
evolution. As demonstrated in Table [5j the matching to this formalism for a variety of dark energy 
and modified gravity models has already been performed. Observational constraints on individual 
models, for which the form of the coefficients is predicted, should be stronger than the general 
bounds that can be set from model-independent analyses. 
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A Applying the Stiickelberg Trick 

In this appendix, we give the transformation law for various quantities using the Stiickelberg trick 
to introduce the it field. 



Functions of time transform as t — > t + it. 



/(*)„2 



f(t) f(t + vr) « f(t) + /(t)TT + + . . . (A.l) 

Such functions will typically be Taylor-expanded about the background. For terms that start at 
second order in perturbations, we ignore the expansion of their time-dependent coefficients, under 
the assumption that they are not rapidly oscillating. 

Operators that are diffeomorphism invariant do not receive any contributions from it. For 
example, the Ricci scalar R. However, when looking at the perturbative terms, e.g. 5R, the 
functions of time forming the background part of the term need to be expanded. For example, 

5R = R- i? (0) -> R - i? (0) - i? (0) 7r - . . . . (A.2) 

Care must be taken when dealing with the normal and extrinsic curvature tensor however, as these 
quantities depend on the foliation, which also changes under a transformation. 

Any tensor that includes explicit time components in the action must transform using the tensor 
transformation law. For example, 

f d{t+ Q x f ]) f" = f° + (A.3) 

We will often need to transform the time-time component of the inverse metric, g 00 (or its pertur- 
bation, which transforms in the same way). It transforms as 



,00 



d(t + 7T(x X ))d(t + 7T(x X )) 00 00 ~ 0i 2 (W) 2 

^ 9 =9 -^ + 2 ^9 +2V^ -vr +— —. 



We will often only need this expansion to first order. 

The normal doesn't just follow the tensor transformation law, as it also depends on the foliation 
of spacetime, which changes with the coordinate transformation. Our foliation is initially described 
by t = const. The normal to these surfaces is given by 

n = M (A.5) 



On the background, this reduces to rffi = 5®. 
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When the tt field is introduced, we have t = t + it, x = x. The index on the normal transforms 
following the tensor transformation law, while the normal itself transforms because of the changing 
foliation. 

^' +T> (A.6, 
Expanded to linear order in tt (and ignoring metric perturbations), this gives 



dx 



-[n^l-^ + d^] , (A.7) 



where n M is the normal to slices of constant time t. 

For the extrinsic curvature tensor, we use the following definition. 

K\ = (g^ + n^)V a n u (A.8) 

Because K^ v n v = K^yn^ = 0, a zero component for the extrinsic curvature tensor will never 
explicitly appear in the unitary gauge action. Thus, the indices for this object will always be 
contracted, and we only need to concern ourselves with the transformation that arises from the 
change to the foliation (and not from the tensor transformation law, which will cancel from the 
transformation of the corresponding contracted index). We use the first order transformation of 
the normal from above. Note that to linear order, no — > uq and n« — > hi + diir, so we should look 
at time and space components separately. 

The extrinsic curvature tensor transforms as follows, to linear order, where we ignore the tensor 
transformation law and drop tildes. 

K% -> K\ + (V A + n CT n A ) V A {-n^ + V^vr) + (V CT ^ - n CT ^) n A V A n M + n CT (v A vr - n A ^) V A n M 

(A.9) 

In terms of individual space and time components, this transforms as 

K\ -+ K\ , (A.10) 
K°i -»• K°i+HViir, (A. 11) 

K^^K^-H^Vjir, (A.12) 
K^K^ + ^VkVjTr. (A. 13) 

CL 

The trace of the extrinsic curvature tensor transforms as 

K%^K% + ^-. (A.14) 



B Metrics and Gauges 

To remain general with respect to the choice of spatial coordinates, we do not specify a spatial 
metric. This is important for working with spatial metrics involving spatial curvature. We write 
the background FRW metric as 

ds 2 = -dt 2 + a(t) 2 g ij dx i dx j , (B.l) 
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where cjij is a three-dimensional spatial metric with no time dependence. To ensure that this 
describes an FRW metric, we require that 



Rijkl = ko {gikQjl - gijQkl) 



(B.2) 



where ko is the curvature constant (we reserve k for wavenumber, and use K for the trace of the 
extrinsic curvature tensor). The metric determinant is 



-9 = a yjg. 



(B.3) 



For metric perturbations, we perform most of our calculations in synchronous gauge, and trans- 
form the results to Newtonian gauge for generality. We consider perturbations of the form 



ds 2 = -dt 2 + a(t) 2 ((jij + hij)dx l dx j . 



(B.4) 



There are two scalar degrees of freedom in hij. The first is defined as the trace h = g^hij = a 2 g lJ hij. 
The second is rj, defined by 



h 



h. 



9ij+ (ViVj-^'VfcVjlij, 



(B.5) 



where V is the covariant derivative associated with g^. As a purely spatial tensor, we will raise 
and lower indices in hij with the spatial metric and its inverse only. 



The connection coefficients can be written to first order in hij as 



1 tt 


= r* ti = r tt = o, 


(B.6) 


r*.. 

n 


= Hgij + Ha 2 hij + -a 2 hij , 


(B.7) 


pi 


= H5) + \ti 3 , 


(B.8) 


-pi 
1 jk 


= f*jfc + \t [v k hji + VjhM - Vihj k . 


(B.9) 



The coefficient r* - fc is written in terms of the spatial metric and its connection r l - fc , which is also 
used by the spatial covariant derivative Vfe. Note that spatial covariant derivatives commute with 
time derivatives. 

Using these connection coefficients, the d'Alembertian of a function of time becomes 



□/ = -/- hH + ^hjf, 

while the d'Alembertian of an arbitrary function to leading order is 



(B.10) 



nf = -f-3Hf + ^V 2 f, 



(B.ll) 



where V 2 represents the spatial Laplacian. 

The curvature invariants to first order in scalar perturbations are given below. We do not 
provide the Riemann tensor here, instead referring readers to Appendix D of Kodama and Sasaki 
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31]. The Ricci tensor is as follows. 

Rtt = 



- 3 (k + H 2 ^j -Hh--h 



Ri 



gij + [H + 3H 2 ) ( gij + a 2 /iy) + -Ha 2 ^ + -a 2 h i3 + -hHg 



2k n 



+ g ( ViVj + g tj V z )[V z 7]-h)+ 2k ViV jr] 



1 • 



Combining these components to calculate the Ricci scalar yields 



R - 6 ( II + 2H 2 + ^ ) + AHh + h - 2^h + I 2k V z r] - ^V z h + -V*V-/> j . 

(XI (X (X \ O O 



kn 



1 



2^2. 



2 ~ 



2~ 



The Einstein tensor becomes as follows. 



G tt 
Gu = Rtt 



3 ( H z + ^ j + Hh- '^h+ I fc V"r/ - -V z /i + -V 2 V 2 r/ 

(X J (X CL \ o o 



I— 



(B.12) 
(B.13) 

(B.14) 

(B.15) 

(B.16) 
(B.17) 



G 



9ij 



3H 2 + 2H+^ 
a z 

1 



a 2 hij ( 3H 2 + 2H + 3^ 



1 



fen 



+ -i/a + -a /ijj - -hHgij - gij-h + ftj^/i 

1, 
6 



+ -V, V ; ( V -77 - /i) - g^2<?iiV z (V^ - h) + 2k ViV jV - h^ 2 V 



(B.18) 

The normal to surfaces of constant time is particularly straightforward in synchronous gauge. 



n. 



yJ-ffPdctdpt " 

We calculate the extrinsic curvature tensor from the normal. We use the definition 

K»p = {51 + n^rf) V a n u . 
To first order in perturbations, the components are 

1 



Knn = Kit 



0, 



Ki 



The trace of the extrinsic curvature K = g^ u K^ u is then given by 

K = -3H - -h. 

2 



Hgij - Ha 2 hij - -a 2 hij. 



(B.19) 



(B.20) 



(B.21) 



(B.22) 



B.l Gauge Transformations 

For some calculations, particularly the anisotropic shear stress and the effective Poisson equation, 
it is convenient to work in Newtonian gauge. The transformations from synchronous to Newtonian 
gauge are well known. See, for example, Kodama and Sasaki [31] or Ma and Bertschinger [22] . 
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Here we provide the formulas for the transformation between synchronous gauge and Newtonian 
gauge. We focus only on the scalar modes of interest. 

Begin with a general metric with perturbations written as 

2 - -(1 + 2^)dt 2 + a{t) 2 [(1 - 2<j>)gij + Xij] dx i dx j , (B.23) 



ds 2 



where Xij is traceless. The special case of ip = defines synchronous gauge, and Xij = defines 



Newtonian gauge. Perform a coordinate transformation x M — > x 
lower indices on e M using the background metric. We choose e M to be 



x^ + e^, where we raise and 



eo 



-a 2 /3. 



(B.24) 



for some scalar function j3. This transformation preserves the form of the perturbed metric. The 
perturbations transform as follows, where we use hats to denote the new quantities. 



{j) = jj, - a 2 (3 - 2a 2 H$ 

= + a 2 Hp + \v 2 p 
1 



Xij = Xij - 2 ( V.Vj - -gijV" ) 



(B.25) 
(B.26) 

(B.27) 



To transform from synchronous to Newtonian gauge, choose (3 = rj/2. The perturbations in 
Newtonian gauge are then as follows. 



N 



1 



a 2 f/ — a 2 Hf] 



2 

h 1 



1 ~ 



+ -cfHri+ -V z rj 

D 



x" = o 



(B.28) 

(B.29) 
(B.30) 



When changing gauge, we must also remember that the n field transforms, as described in Appendix 
|a| Under the transformation described above, we have tt = ir — a 2 $, which becomes 



7r N = n s 



-V, 



(B.31) 



when transforming from synchronous gauge to Newtonian gauge. 

When transforming tensors, the tensor transformation law also applies. For example, the Ein- 
stein tensor transforms as the following to linear order. 



= C + G a e M 



G^ e a 



(B.32) 



In momentum space, the perturbations in synchronous and Newtonian gauge are related by 

,2 



V 



N 



k 2 



1 

2 

a 2 H 



(h + 677) + H(h + 677) 
h + 6fj 



2k 2 



7T 



7T + 



2k 2 



(h + 6^) . 



(B.33) 
(B.34) 
(B.35) 



Note that these equations use f\ rather than r\. It is also useful to note that ir s = 7T — ip N . 
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C Momentum Space 



In this appendix, we detail the formalism we will use for the transformation to momentum space. 
This formalism is necessary for working with a metric involving spatial curvature. 

We follow the conventions of Kodama and Sasaki [31j . For each k, define Y^(x l ) to be a solution 
of the equation 

-k 2 Y r , 



V 2 Y, 



(C.l) 



From now onwards, we suppress the k dependence of Y. Both Y and Y* will be solutions of this 
equation, as will any linear combination. We fix this freedom by choosing 



lim Y{x % 
fc -K) 



Ak-x 



(C.2) 



so that the modes become the usual Fourier modes when the background spatial metric is flat. We 
choose the normalization of Y such that 

(C.3) 

(C.4) 
(C.5) 



tPxy/gYffi = (27T) 3 5 3 (k-k'). 

Taking derivatives of Y, we define vector and tensor mode functions as 

Yi = -k- 1 \I l Y , 
Yij = k- 2 ViVjY + -g i:j Y. 



Like the metric perturbation hij, we raise and lower indices on Y{ and Yy with the spatial metric 
only. This is useful not only for consistency with hij, but also because it ensures that Y, Yi, and 
Y^ are independent of time, no matter the position of their indices. 

In synchronous gauge, the metric perturbation decomposes as 



hij {x, £) 



d A k 



h(k,t) 



Ygij + k 2 ri(k,t)Yi 



Following Ma and Bertschinger [22], we define a new field rj in momentum space by 

1 



i](k,t) 



k 2 



[h(k,t) + 6fj(k,t)). 



Using this definition, the metric perturbation in synchronous gauge becomes 



hi 



d s k 



^Y~ 9lJ -(h + Gf})Yij 



(C.6) 



(C.7) 



(C.8) 



In the limit of no spatial curvature, the mode function Y = exp(ik-x) yields Ma and Bertschinger's 
Eq. (4). 



D Equations of Motion 

Here, we present the contribution to the n equation of motion and the effective stress-energy tensor 
for a number of operators, both in synchronous gauge and Newtonian gauge. We begin by giving 
the equations for the background operators, then present the contribution each further operator 
adds to these equations. 
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D.l Background Terms 
D.l.l 7r Equation of Motion 

In synchronous gauge, the equation of motion for ir from the background action [first line of Eq. 



(3.9)] is 



= ctt + c^Jvr + + 3Hck + ( - 3Hc) vr + C H - ^QSR. 
a 2 \ 4 J 2 4 

The perturbation to the Ricci scalar, included in the final term, is given by 

5R = 4Hh + h + \f) (3fc - k 2 ) . 



(D.l) 



(D.2) 



In the limit of c(t) = 0, the equation of motion for tt becomes algebraic in the absence of further 
operators, and can be solved to give 

5R 



IT 



(D.3) 



However, the {ir, h, ff\ system still contains one dynamical degree of freedom because of the kinetic 
mixing. 

In Newtonian gauge, the equation of motion from the background action is given by 



= cn N - c^ N + £n N + c(n N - i> N ) + 3cHir N - 3c<j) - 6cH^ N 



-V 



mr 



nsR 



N 



(D.4) 



where quantities in Newtonian gauge are denoted by a superscript N. In Newtonian gauge, SR N 
is given by 



6R 



N 



a 2 

N ■ 



k 2 

i.N i n K J.N 



N 



1N A (k 2 - 3k ) 



■N 



(D.5) 



In the appropriate limit, Eq. (D.3) still holds, using 5R instead of 5R. 



D.l. 2 Effective Stress-Energy Tensor 



The gravitational equations of motion can be obtained by varying the action (3.3) with respect 



to the metric. The resulting equations can be written in terms of the Einstein tensor, matter 
stress-energy tensor, and an effective dark energy stress-energy tensor as 



(D.6) 



We ignore the background contributions that have already been solved. The perturbed Einstein 
tensor is given by the following in synchronous gauge, where we use mode functions as described 
in Appendix [O 



6G% 
5G\ 



-Hh + {2k 2 - 6fc ) 
2rjk-(h + 6rj)~ Y 



Y 



k 

Ltt 2k2 ~ 
hH -\ =n 

3 3a 2 1 



h 



2k „ 



S > Y+ 2 



-3H(h + 6fj) -{h + 6fj) + 



2k 2 



Y\ 



(D.7) 
(D.8) 

(D.9) 
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The scalar modes of the matter stress-energy tensor in momentum space can be decomposed in the 
following manner, following Kodama and Sasaki [31J. 

(D.10) 
(D.ll) 
(D.12) 



T° {m) = -(Pm + SpY) 
T\ {m) = {p m + P m )vY l 
T) {m) = (P m + SPY) 5) + P m nY 



Working to linear order, the components of Eq. ( D.6[ ) are as follows, where we again use synchronous 
gauge and focus on terms arising from the background operators only. 
Time-time (Synchronous) 



2mf 



n 



k f] 



a 2 H 



-h — 3kofj 



5p - pqtt - 2ck 



(D.13) 



Time-space (Synchronous) 



2k 2 fj - k (h + 6fj) = {p m + P m )kv + mfotf-k + [p Q + P Q ] k 2 n 



(D.14) 



Space-space trace (Synchronous) 



h ■ 2k 2 2k 



5P + Pqtt + (p Q + P Q ) tt 



n 



2 k 2 h 



n + ii + 3-Htt + -^vr + - + vr ( 3H Z + 
3 a z 3 



(D.15) 



Space-space traceless (Synchronous) 

~2k 2 



-fj - 3H(h + 6fj) -(h + 6fj) 



P m U + mltl 



7T + — ( h + 677 



(D.16) 



We can translate these into Newtonian gauge. The required gauge transformations are detailed 
in Appendix [B] 

Time-time (Newtonian) 



2k 2 - 6k 



N 



+ 6H(</> N + Hip 



,,N\ 



5p N -p Q ^ -2c(t)(n» -r) 



JV 



N „i.N\ 



_.\ / rj2 r'l 1 . nrjl-N j.N\ 1 k 2 N 



3k n [H 2 - H + ^ + 3H (tt n - V w ) + 



~ ■ 3(o v • II r x 



(D.17) 



Time-space (Newtonian) 



mln(t)2k 2 



(p m + P m )kv N + [p Q + Pq] k 2 n N + m 2 nk 2 (ir N - V W ) (D.18) 
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Space-space trace (Newtonian) 



2k 2 

2<P + 2 ( 3H 2 + 2H ) if, + 2H(iP + 30) + ^ (0 " 



^ v +PQ7r JV + (pQ + P Q )(^ 

"ft 
f2 



»/.JV 



20^ + 3^ 



5^ + ^(3^ + ^ 



+ 



2 k 2 
3a 2 



f7T 



Space-space traceless (Newtonian) 



k 2 



a* a 
The matter terms also transform in going to Newtonian gauge 



5p N = 5p + p m (h + 6fj) 



2k 2 



5P n =5P + p m _(/ l + 67?) 



(D.19) 
(D.20) 

(D.21) 
(D.22) 
(D.23) 



(p m + P m )kv = (p m + P m )kv + (P m + p m ) y + 
The matter anisotropic shear stress is invariant. 

D.2 Further Operators 

We now investigate the terms (5g 00 ) 2 , 5g 00 5K 1 "^, (SK^) 2 , 5K^ U 5K U ^ and 5g 00 5R^\ calculating 
their contributions to the equations of motion. 



To describe contributions to the tt equation of motion (D.l), we write the equation of motion 



as 



= CTT + C^j7T + . . . + A,,- , 



(D.24) 



and state A„ for individual operators, where the ellipsis includes terms from other operators. For 
the stress-energy tensor, we state the contribution to T^ v , for which there should be no ambiguity. 



Operator: (5g 00 ) 2 

This term very generically arises from writing a covariant theory in our formalism. As a simple 
example, any model that involves f(X) type terms where X oc g^ v V u(f>V v <j) such as fc-essence will 
include this term. After implementing the Stiickelberg trick, the action becomes the following. 



5, 



(^00)2 



(Sg 00 - 2*) 1 



The 7r equation of motion receives the following contribution in synchronous gauge. 

Ajf = 2M 2 4 # + (2M 2 4 + 6HM$)-k 



(D.25) 



(D.26) 
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In Newtonian gauge, this is 

= 2M^ N - i) N ) + (2M 2 4 + 6HM*)(ir N - ij> N ). (D.27) 
The contribution to the stress-energy tensor is 

K (Q W>o )2 = -^Ml^SlY = -4M 2 4 (^ - (D.28) 
in synchronous gauge and Newtonian gauge respectively. 

Operator: Sg^SK^ 

This term also commonly arises from writing covariant theories in this formalism. For example, 
kinetic braiding and galileon terms will generate it, while the full expansion of quadratic curvature 
invariants in this formalism will also typically generate it. 

The action with the tt field is 



S 5g 00 5K = I d X 



(D.29) 



7r equation of motion in synchronous gauge: 



1 h ■ k 2 7r\ M? ( h ■■ k 2 ir 



In Newtonian gauge: 

A w = 1 -[{W + d t )Mf] (3 (> + H^ N ) + 3H*» - ^tt n 

+ ^1(34," + 6H^ N + 3H^ N + 3Hir N + 2H~k n - ~i> N ) (D.31 i 

2 \ a z a z 

Stress-energy tensor (synchronous gauge): 

T °o iQ) 6 9 oo SK = 'Ml (~\ + 321* - ^vr - 3H*j Y (D.32) 

T V Q) 8 9 »o SK =kM?nY t (D.33) 
T ) {Q) 5 g oo SK = (322 + d t ) [Mfn] b)Y (D.34) 

Note that we find the ij component is a factor of two different to Eq. (C.13) in Creminelli et al 
[6]. In Newtonian gauge: 



< iQ Lo OSK = -Mf hU N + H^) + 3H* N -~ 2 n N -3H U N — tj) N ) ) Y (D.35) 



A: 2 

T^ Q \ g005K =kMl{^-^)Y i (D.36) 
T ) {Q \ g oo 5K = (3H + ft) (tt* - ^)] <^ (D.37) 
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Operator: (SK) 2 



The effect of this term on the behavior of the tt field has been investigated in detail in the ghost 
condensate limit, in particular in Ref. [3J. 

Action with tt field: 



\8K) 2 



d xy/—g 



Mf(i) 



K^„ + 3H + 3Htt + 



V 2 vr 



tt equation of motion in synchronous gauge. 



3Htt 



.-7T 



In Newtonian gauge: 



2 ozj * \ o IN , ttj.N\ . oiV_A r K _AT 



2 V a 
Stress-energy tensor in synchronous gauge: 



+ Hifr + 3HTT 1 
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— + 3Htt kTT 



8}Y 



(D.38) 



(D.39) 



(D.40) 



(D.41) 
(D.42) 
(D.43) 



Again, we find the ij component to be a factor of two different from Creminelli's Eq. (C.13). 
Newtonian gauge: 



^0 (Q) 
{SKf 



-:\IIM 2 ( 3 (j> N + Hi, N ) + 3^ - ^ 



TT N \Y 



,o (Q) 



-kM% ( 3 ( <j> N + Hip N ) + 3#V 
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+ 3Htt 
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;TT 
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(D.44) 
(D.45) 
(D.46) 



Operator: 5K^ v 8K h ^ 

In terms of the effect on the behavior of the tt field, this term is qualitatively similar to (5K ) 2 . The 
effect on the gravitational degrees of freedom is rather different, however. Action with tt field: 



d^x^J^g 



2 

2 



~ik 



5K) + Htt5) + ^VfcVjTT ) ( 5K\ + HnSj + ^ V/V.vr 



(SK\y + 2 (6K\ - H^VjTt ) (SlCi + HViTT 
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(D.47) 



The 7r equation of motion in synchronous gauge is 

k k , TT o • \ 3 



A. = M 2 



2 i 



( kT fe" /TT o 3 - 9 \ / k 2 H\; k , 



and in Newtonian gauge: 



A„ = M 2 



( k 4 k 2 , rr9 , tVn , 3 tV2 \ n _ { k 2 3H 



V2a 4 a 



2 (H 2 + H) + -H')n I 



\2a 2 2 / 



The stress-energy tensor in synchronous gauge is 



T (Q) 

J SK^SK"^ 
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T (Q) 

1 i SK^SK"^ 


= MlkYi 


rpi (Q) 
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(D.49) 

(D.50) 
(D.51) 



Ma / J J 



' A: 2 



or a 

2 ; 2 



l riv f* 2 . A; 2 / 2M 3 \ 1.. /3 rr M 3 



a 



- ' ' ' " 1 :>/;: y " 2" ! V2" : /•..) / '} 



M 3 , 
2 " f M 3 - 



and in Newtonian gauge, 
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(D.52) 

(D.53) 
(D.54) 

.N 



A: 2 



,1 

2Ms\ 
M 3 



(D.55) 



Operator: Sg 00 8R^ 

This term is required to fully describe perturbations in Horndeski's theory. Note that is the 
Ricci scalar of the spatial metric g^ , without a prefixed scale factor a 2 . 

The action with the ir field introduced is 



S Sg 00 SR( 3 ) - J 



M ^> (Sg 00 -2n) 6RW 



The perturbed three-dimensional Ricci scalar is 

5RW = (h + 6fj)2k - 4k 2 fj = -4k 2 <f> N 
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(D.56) 



(D.57) 



in synchronous and Newtonian gauge, respectively. 

The 7r equation of motion is as follows, in synchronous and Newtonian gauge respectively. 



1 



1 



A n = - [3HM 2 + 2MM) [(h + 6fj)2k - Ak 2 fj] + -M 2 (t) {h + 6^)2£;o - ^k 2 fj 



-2 [3HM 2 + 2MM) k 2 cj> N - 2M 2 (t)k 2 <j) N 



The stress-energy tensor in synchronous gauge is 
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In Newtonian gauge, this becomes 
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